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1. INTRODUCTION 

The discrepancies between the values exactly obtained on the variables under consideration for sampled units 
and the corresponding true values are termed as measurement errors. In general, standard theory of survey 
sampling assumes that data collected through surveys are often assumed to be free of measurement or response 
errors. In reality such a supposition does not hold true and the data may be contaminated with measurement 
errors due to various reasons; see, e.g., Cochran (1963) and Sukhatme et cil (1984). 

One of the major sources of measurement errors in survey is the nature of variables. This may happen in case of 
qualitative variables. Simple examples of such variables are intelligence, preference, specific abilities, utility, 
aggressiveness, tastes, etc. In many sample surveys it is recognized that errors of measurement can also arise 
from the person being interviewed, from the interviewer, from the supervisor or leader of a team of interviewers, 
and from the processor who transmits the information from the recorded interview on to the punched cards or 
tapes that will be analyzed, for instance, see Cochran (1968). Another source of measurement error is when the 
variable is conceptually well defined but observations can be obtained on some closely related substitutes 
termed as proxies or surrogates. Such a situation is encountered when one needs to measure the economic status 
or the level of education of individuals, see Salabh (1997) and Sud and Srivastava (2000). In presence of 
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measurement errors, inferences may be misleading, see Biemer et al (1991), Fuller (1995) and Manisha and 
Singh (2001). 

There is today a great deal of research on measurement errors in surveys. An attempt has been made to study 
the impact of measurement errors on a family of estimators of population mean using multiauxiliary 
information. 

2. THE SUGGESTED FAMILY OF ESTIMATORS 

Let Y be the study variate and its population mean |i (l to be estimated using information on p(>l) auxiliary 
variates Xi, X 2 , ...,X p . Further, let the population mean row vector /d ' — (jd l , /I 7 , • • • , ) of the vector 

x' = (x u x 2 ,x p ). Assume that a simple random sample of size n is drawn from a population, on the study 

character Y and auxiliary characters Xi, X 2 , ...,X P . For the sake of simplicity we assume that the population is 
infinite. The recorded fallible measurements are given by 



yj= y j +E j 

X ii = X;; 1 = 1 , 2 , 1 

j = 1,2, •••,«. 

where Yj and Xy are correct values of the characteristics Y and X, (i=l,2,..., p; j=l,2,..., n). 

For the sake of simplicity in exposition, we assume that the error Ej's are stochastic with mean 'zero' and 
variance a (0 ) 2 and uncorrelated with Yj's. The errors r|ij in xy are distributed independently of each other and of 
the Xy with mean 'zero' and variance o (i) 2 (i=l,2,...,p). Also Ej's and T|y's are uncorrelated although Yj's and Xy's 
are correlated. 

Define 



M, =— , (i =1,2 

Hi 



u 



(u 1 ,u 2 ,---u p ) Up ,e T =(l,l, •••,!) 



1 Xp 



1 n 

y-zly, 



1 n 

x,. =-^x. 



n 



ij 



j = t 



With this background we suggest a family of estimators of p 0 as 

H g =g(y,u T ) 
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( 2 . 1 ) 



where g ( y, U T ) is a function of y,U 1 ,U 2 ,"-,U p such that 

8(u 0 S) = V o 

_ dg(-) 



dy 



= 1 



(/V 



and such that it satisfies the following conditions: 

1. The function g(y,U T ^ is continuous and bounded in Q. 

2. The first and second order partial derivatives of the function g(y,u T ) exist and are continuous and bounded 
inQ. 

To obtain the mean squared error of /} „ , we expand the function g{y,U T ^j about the point (p 0 e T ) in a second 
order Taylor's series. We get 



v g = g(vo’ eT )+(y-v 0 f^ 






+ (u-e) T g% 0 ,e T ) 



+ — 
2 






+ 2 (y - Vo \ u - e )‘ 



3g (1) Q 



dy 



+ 



(u - e)' g <2) (y*, u * T \u - <?)) 



(2.2) 



where 



y* = Ho+d(y- fl Q \ u* = e + e(u - e), 0 < 9 < 1; g (1, (-) 

denote the p element column vector of first partial derivatives of g(-) and g <2> (-) denotes a pxp matrix of second 
partial derivatives of g(-) with respect to u. 

Noting that g(p 0 .e T )= Po. it can be shown that 

E{fi g )=ii 0 +C>(«- 1 ) 

(2.3) 

which follows that the bias of fL g is of the order of n" 1 , and hence its contribution to the mean squared error of 
/} will be of the order of n" 2 . 



From (2.2), we have to terms of order n" 
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MSE (fl g )= E{y - p 0 )+ (u - e) T g (1) (, 0 ,/)} 2 

= E i{y = p 0 f + 2 (y - p 0 X» - e Y 8 (1 W) 

+ (g % 0 -e T )) T (« - e\u - ef (g ( % 0 , e r )). 

= -kfc + C ( o } )+ 2fl 0 b T g% o , e T )+ (g { %S)J A(g % 0 S)) 

n 



where b T =(bi,b 2 ,...,b p ), bi,=p 0i C 0 Ci,(i=l,2, 
Ci=Oi/|ii, C (i )= Oi/p„ (i=l,2, ...,p) and C 0 =Oo/|Xo, 





cf + q 2 ) 


Pl2 E l E 2 


Pi 3 CtC 3 


•• P lp c l c p 




p n C l C 2 


C 2 + C ( 2 ) 


p 23 C 2 C 3 


P2 p ClC p 


A = 


p l3 C x C 2 


p 23 C 2 C 3 


c 3 2 + c ( 2 3) • 


•• P2 P C,C P 




_P lp C 1 C p 


P2„C 2 C p 


Pi p C 3 C p • 


•• c >c 2 p) 



The mse(uJ at (2.4) is minimized for 

g { %S) = -H 0 A~ 1 b 



Thus the resulting minimum MSE of p ? is given by 

min.MSE (fi g )= (u 2 / n\cl + C 2 0) -b r A~ l b\ 



Now we have established the following theorem. 
Theorem 2.1 = Up to terms of order n , 



(2.4) 



(2.5) 



( 2 . 6 ) 



MSE(f u g )> (p 2 !n\cl +C{ 0) -b T A- l b\ 



(2.7) 



with equality holding if 



g l \s) = -fi 0 A- 1 b 



It is to be mentioned that the family of estimators f l g at (2.1) is very large. The following estimators: 

p 

I>, = 1 ’ 



K {1) = yl<o 






v x '-y 



[Olkin (1958)] 
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( 2 ) _ TT 



/V ' = 






vZ®- 

1=1 

P 

X®^' 



' X ' 



K pi J 



r 

X £ »/ [Singh (1967)] 



y^- 

^ n 



/>. w = 



x®* 

1=1 

p 

x®.*. 



-, ^CJ ( . = 1, [Shukla (1966) and John (1969)] 



v-^ 

^ n 



X®-^ 



-; 2®- =!’ [Sahai et al (1980)] 



1=1 



= , 2^,^/ = I 5 [Mohanty and Pattanaik (1984)] 

i = 1 X i 1=1 



/ \®i 

Pi 



y x , j 



p* {6) = y 



/i ? (7) = yPJ— , 2« ; =1, [TntejaandSahK!"!)] 
1=1 1=1 



'f °h x , ' 

M Pi 



I' 

, 2®,- = 1) [Mohanty and Pattanaik (1984)] 

1=1 



JL 



~ ( 8 ) - 

/v = y 



£, (9) = y 



/V ( 10 ) _ 

/V -y 



/V (11) _ 

/V = y 



,= i y 



^ m iPi 

i = 1 X i 



®^i + X®.- 

f=l 



, = 1) [Tuteja and Bahl (1991)] 

p+i 

> X®I =1 - 



1=1 






v x ‘ yj 



® P+ i+X®t 






yj 



i=l 



p+i 



, X®t =L 
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X®/ + X®.- 


l X i J 


i=q+l 


J- 




*=1 1=<7+1 ^ 



=1 



; [Srivastava (1965) and Rao 
and Mudhalkar (1967)] 



fL „^ 12) = yTT — (a, ’s are suitably constants) [Srivastava (1967)] 

6 X A , , 



r \ a > 



i=i \P, ; 



P 



( 13 ) _ T7 



yfl 2 - 



r~ \ a ‘ 

x. 



v^y 



[Sahai and Rey (1980)] 
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M k + <*,(*, -Mi)} 



[Walsh (1970)] 



/i l , (l5) = y exp jy 1 0. log M, | [Srivastava (1971)] 



/},, (l6) = yexpjy 1 — l)| [Srivastava (1971)] 



- y^j m i exp{(0. / (Of )logW,.}, ^j(Of = 1, [Srivastava (1971)] 



/*, p8) = y+X«,ft-A 0 



etc. may be identified as particular members of the suggested family of estimators fl . The MSE of these 
estimators can be obtained from (2.4). 



It is well known that 



V(j) = (i^o / n \Co + C ( 2 0) ) 



( 2 . 8 ) 



It follows from (2.6) and (2.8) that the minimum variance of is no longer than conventional unbiased 



estimator y . 

On substituting O(0) 2 =0, a ( i) =0 Vi=l,2,...,p in the equation (2.4), we obtain the no-measurement error case. In 
that case, the MSE of , is given by 

MSE(u ? )= -[clill + 2 ^b T g *% 0 j) +(g A*{g * { % 0 ,/)) 

n 

= MSE (fi *) 



(2.9) 



where 



V* =8 



= (7 * 



*9 9 9 9 



V 



A*t li 2 H, 



= g*iy,u T ) 



J 



(2.10) 
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and Y and Xj (/ = 1,2, •••,/?) are the sample means of the characteristics Y and X, based on true 
measurements. (Yj,Xjj, i=l,2,...,p; j=l,2,...,n). The family of estimators fl, A * at (2.10) is a generalized 
version of Srivastava (1971, 80). 

The MSE of fL g * is minimized for 

g*% 0 ,eT) = -A*- 1 b/l 0 

(2.11) 

Thus the resulting minimum MSE of fl * is given by 



min.MSE(/i ? *) 




»] 




n 



('-«=) 



( 2 . 12 ) 

where A*=[a*;j] be a pxp matrix with a*^ = pijQCj and R stands for the multiple correlation coefficient of Y on 

From (2.6) and (2.12) the increase in minimum MSE (/ i ? ) due to measurement errors is 
obtained as 



min 



.MSE(/2 )-min.MSE(/2 *) = 



f 



^\c( 0) + b T A*- l b-b T A- l b\ 



V 

>0 



This is due to the fact that the measurement errors introduce the variances fallible measurements of study variate 
Y and auxiliary variates Xj. Hence there is a need to take the contribution of measurement errors into account. 



3. BIASES AND MEAN SQUARE ERRORS OF SOME PARTICULAR ESTIMATORS IN THE 
PRESENCE OF MEASUREMENT ERRORS. 

To obtain the bias of the estimator /} , , we further assume that the third partial derivatives of g(y, U T ) also 

exist and are continuous and bounded. Then expanding g(y,li T ) about the point (y , li T ) = (jLL 0 , e T ) in a 
third-order Taylor's series we obtain 






= g(Mo’ eT ) + (y-Vo} 



M) 

dy 



l(/V' 



+ {u-e) T g% 0 ,e T ) 
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+2(v-/a»-*)vw) 

+ (u-e) T (g { % 0 ,e T )\u-e)} 

+ ^{(y-^°)^+(M-^)^-} g(f>u* T ) 

(3.1) 

where g‘ l2l (|Io,c T ) denotes the matrix of second partial derivatives of g (y, li T ) at the point 

(y,u T )=(ju 0 ,e T ). 

Noting that 

g(u 0 e T )=/u 0 

Ml =i 




values of g (1) (|i 0 ,e T ), g f2l C|_L 0 ,e T ) and g <l2l (|Io,e T ) are given in the Table 3.1. 




Table 3.1 Biases and mean squared errors of various estimators of |i« 
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Table 3.1 Biases and mean squared errors of various estimators of |i« 




where OC * T =(oc 1 *,a 2 *,...,cXp*) with OL * ( =(cx i ,(a 1 , i=l,2,...,p) 
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4. ESTIMATORS BASED ON ESTIMATED OPTIMUM 

It may be noted that the minimum MSE (2.6) is obtained only when the optimum values of constants involved in 
the estimator, which are functions of the unknown population parameters p 0 > b and A, are known quite 
accurately. 

To use such estimators in practice, one has to use some guessed values of the parameters |i 0 , b and A, either 
through past experience or through a pilot sample survey. Das and Tripathi (1978, sec. 3) have illustrated that 
even if the values of the parameters used in the estimator are not exactly equal to their optimum values as given 
by (2.5) but are close enough, the resulting estimator will be better than the conventional unbiased estimator y . 
For further discussion on this issue, the reader is referred to Murthy (1967), Reddy (1973), Srivenkataramana 
and Tracy (1984) and Sahai and Sahai (1985). 

On the other hand if the experimenter is unable to guess the values of population parameters due to lack of 
experience, it is advisable to replace the unknown population parameters by their consistent estimators. Let (p 
be a consistent estimator of (f)=A‘ 1 b. We then replace cf> by <p and also p 0 by y if necessary, in the optimum 
g resulting in the estimator fi g ^ est j , say, which will now be a function of y , u and (f>. Thus we define a 
family of estimators (based on estimated optimum values) of p 0 as 



Q ^ ^ 



Hg(est) 8 



(y,“ T J T ) 



where g**(-) is a function of 



(y,u T ,(p T ) 



such that 



g**(u 0 ,e T y)= for all fi 0 , 

d g **(.)! 



dy 



= 1 



{fl 0 ,e T ,<l> T 



(4.1) 



dg **0 



du 



3*0 






du 



= -Ho A ~' b = 



(Mo ,ef 



(4.2) 



and 



dg**Q 

dip 



(Ho.e 4 



= 0 
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With these conditions and following Srivastava and Jhajj (1983), it can be shown to the first degree of 
approximation that 



MSE^ M )=minMSE(/i g ) 
[ Ct + C? 0) -b T A-'b\ 






V J 

Thus if the optimum values of constants involved in the estimator are replaced by their consistent estimators and 
conditions (4.2) hold true, the resulting estimator fi g ^ esl ^ will have the same asymptotic mean square error, as 



that of optimum . Our work needs to be extended and future research will explore the computational aspects 
of the proposed algorithm. 



REFERENCES 



Biermer, P.P., Groves, R.M., Lyberg, L.E., Mathiowetz, N.A. and Sudman, S. (1991): Measurement errors in 
surveys, Wiley, New York. 

Cochran, W. G. ( 1963): Sampling Techniques, John Wiley. New York. 

Cochran, W.G. (1968): Errors of measurement in statistics. Technometrics, 10(4), 637-666. 

Das, A.K. and Tripathi, T.P. (1978): Use of auxiliary information in estimating the finite population variance. 
Sankhya, C, 40, 139-148. 

Fuller, W.A. (1995): Estimation in the presence of measurement error. International Statistical Review, 63, 2, 
121-147. 

John. S. (1969): On multivariate ratio and product estimators. Biometrika, 533-536. 

Manisha and Singh, R.K. (2001): An estimation of population mean in the presence of measurement errors. 

Jour. Ind. Soc. Agri. Statist., 54 (1), 13-18. 

Mohanty, S. and Pattanaik, L.M. (1984): Alternative multivariate ratio estimators using geometric and harmonic 
means. Jour. Ind. Soc. Agri. Statist., 36, 1 10-118. 

Murthy, M.N. (1967): Sampling Theory and Methods, Statistical Publishing Society, Calcutta. 

Olkin, I. (1958): Multivariate ratio estimation for finite population. Biometrika, 45, 154-165. 

Rao, P.S.R.S. and Mudholkar, G.S. (1967): Generalized multivariate estimators for the mean of a finite 
population. Jour. Arner. Statist. Assoc. 62, 1009-1012. 

Reddy, V.N. and Rao, T.J. (1977): Modified PPS method of estimation, Sankhya, C, 39, 185-197. 

Reddy, V.N. (1973): On ratio and product methods of estimation. Sankhya, B, 35, 307-316. 

15 




Salabh (1997): Ratio method of estimation in the presence of measurement error. Jour. Ind. Soc. Agri. Statist., 



52, 150-155. 

Sahai, A. and Ray, S.K. (1980): An efficient estimator using auxiliary information. Metrika, 27, 271-275. 

Sahai, A., Chander, R. and Mathur, A.K. (1980): An alternative multivariate product estimator. Jour. Ind. Soc. 
Agril. Statist., 32, 2, 6-12. 

Sahai, A. and Sahai, A. (1985): On efficient use of auxiliary information. Jour. Statist. Plann. Inference, 12, 203- 

212 . 

Shukla, G. K. (1966): An alternative multivariate ratio estimate for finite population. Cal. Statist. Assoc. Bull., 

15, 127-134. 

Singh, M. P. (1967): Multivariate product method of estimation for finite population. Jour. Ind. Soc. Agri. 
Statist., 19, (2) 1-10. 

Srivastava, S.K. (1965): An estimator of the mean of a finite population using several auxiliary characters. Jour. 
Ind. Statist. Assoc., 3, 189-194. 

Srivastava, S.K. (1967): An estimator using auxiliary information in sample surveys. Cal. Statist. Assoc. Bull., 

16, 121-132. 

Srivastava, S.K. (1971): A generalized estimator for the mean of a finite population using multiauxiliary 
information. Jour. Amer. Statist. Assoc. 66, 404-407. 

Srivastava, S.K. (1980): A class of estimators using auxiliary information in sample surveys. Canad. Jour. 
Statist., 8, 253-254. 

Srivastava, S.K. and Jhajj, H.S. (1983): A class of estimators of the population mean using multi-auxiliary 
information. Cal. Statist. Assoc. Bull., 32, 47-56. 

Srivenkataramana, T. and Tracy, D.S. (1984):: Positive and negative valued auxiliary variates in Surveys. 
Metron, xxx(3-4), 3-13. 

Sud, U.C. and Srivastava, S.K. (2000): Estimation of population mean in repeat surveys in the presence of 
measurement errors. Jour. Ind. Soc. Ag. Statist., 53 (2), 125-133. 

Sukhatme, P.V., Sukhatme, B.V., Sukhatme, S. and Ashok, C. (1984): Sampling theory of surveys with 
applications. Iowa State University Press, USA. 

Tuteja, R.K. and Bahl, Shashi (1991): Multivariate product estimators. Cal. Statist. Assoc. Bull., 42, 109-1 15. 
Tankou, V. and Dharmadlikari, S. (1989): Improvement of ratio-type estimators. Biom. Jour. 31 (7), 795-802. 
Walsh, J.E. (1970): Generalization of ratio estimate for population total. Sankhya, A, 32, 99-106. 



16 




